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ABSTRACT 

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18deg 2 medium- 
deep survey at 3.6 and 4.5 /J.m with the post-cryogenic Spitzer Space Telescope to ~ 2 /xJy (AB = 23.1) 
depth of five highly observed astronomical fields (ELAlS-Nl, ELAIS-S1, Lockman Hole, Chandra 
Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a 
function of environment from z ~ 5 to the present day, and is the first extragalactic survey that is 
both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters 
at z <; 1 into their cosmological context. SERVS is designed to overlap with several key surveys at 
optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented 
view of the formation and evolution of massive galaxies. In this article, we discuss the SERVS survey 
design, the data processing flow from image reduction and mosaicking to catalogs, and coverage of 
ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science 
results from the survey. 

Subject headings: Astrophysical data, surveys 



1 Infrared Processing and Analysis Centcr/Spitzer Science Cen- 
ter, California Institute of Technology, Mail Code 220-6, Pasadena, 
CA 91125, USA 

2 National Radio Astronomy Observatory, 520 Edgcmont Road, 
Charlottesville, VA 22903, USA 

3 Department of Physics and Astronomy, University of Sussex, 
Falmer, Brighton, BN1 9QH, UK 

4 Center for Astrophysics Research, University of Hertfordshire, 
Hatfield, ALIO 9AB, UK 

5 Institute of Cosmology and Gravitation, University of 
Portsmouth, Dennis Sciama Building, Burnaby Road, Portsmouth, 
POl 3FX, UK 

6 Department of Astronomy, Univcrsita di Padova, Vicolo 
deH'Osscrvatorio 3, 35122, Padova, Italy 

7 Astrophysics Group, Physics Department, University of the 
Western Cape, Private Bag X17, 7535, Bellville, Cape Town, South 
Africa 

8 Department of Physics, University of California, One Shields 
Avenue, Davis, CA95616, USA 

9 Institute of Astronomy, University of Cambridge, Madinglcy 
Road, Cambridge, CB3 0HA, UK 

10 National Optical Astronomy Observatory, 950 North Cherry 
Avenue, Tuscon, AZ 85719, USA 

11 Observatorio Astronomico de Lisboa, Faculdade de Cicncias, 
Univcrsidadc de Lisboa, Tapada da Ajuda, 1349-018 Lisbon, Por- 
tugal 

Centro de Astronomia e Astroffsica da Universidade de Lis- 
boa, Tapada da Ajuda, 1349-018 Lisbon, Portugal 

13 Department of Physics and Astronomy, University of 
California-Riverside, 900 University Avenue, Riverside, CA 92521, 
USA 

14 Department of Physics, McGill University, Ernest Rutherford 
Building, 3600 rue University, Montreal, Quebec H3A 2T8, Canada 

15 International Centre for Radio Astronomy Research, Univer- 
sity of Western Australia, M468, 35 Stirling Highway, Crawley WA 
6009, Australia 



16 Commonwealth Scientific and Industrial Research Organisa- 
tion, Astronomy & Space Science, PO Box 76, Epping, NSW, 1710, 
Australia 

17 Department of Physics, Drexel University, 3141 Chcsnut 
Street, Philadelphia, PA 19014, USA 

18 Institute of Geophysics and Planetary Physics, Lawrence Liv- 
ermore National Laboratory, 7000 East Ave., Livermore, CA94550, 
USA 

19 Department of Physics, University of Durham, South Road, 
Durham, DH1 3LE, UK 

20 Institute for Astronomy, University of Edinburgh, Royal Ob- 
servatory, Blackford Hill, Edinburgh, EH9 3HJ, UK 

21 Institutio de Astroffsica de Canarias, C/Vi'a Lactca s/n, 
38200, La Laguna, Tenerife, Spain 

22 Oxford Astrophysics, Denys Wilkinson Building, Kcble Road, 
Oxford, OX1 3RH, UK 

23 Astrophysics Group, Blackett Laboratory, Imperial College, 
Prince Consort Road, London, SW7 2BW, UK 

24 Dipartimento di Scienze Fisiche, Universita Federico II and Is- 
tituto Nazionalc di Fisica Nuclcare, Sczione di Napoli, Complesso 
Universitario di Monte S. Angclo, Via Cintia, Edificio 6, 1-80126 
Napoli, Italy 

25 Institute for Computational Cosmology, Durham University, 
South Road, Durham, DH1 3LE, UK 

26 Space Telescope Science Institute, 3700 San Martin Drive, 
Baltimore, MD 21218, USA 

27 Department of Earth Sciences, National Taiwan Normal Uni- 
versity, N°88, Tingzhou Road, Sec. 4, Taipei 11677, Taiwan 
(R.O.C.) 

28 Institute for Astronomy, University of Hawaii, 2680 Wood- 
lawn Drive, Honolulu, HI 96822, USA 

29 Institute for Space Imaging Science, University of Calgary, AB 
T2N 1N4, Canada 

30 European Southern Observatory, Karl-Schwartzschild-Strassc 
2, 85748, Garching, Germany 



2 



Mauduit et al. 



1. INTRODUCTION 



Progress in extragalactic astronomy has been greatly 
enhanced by deep surveys such as the Great Observa- 
tories Origins Deep Survey (GOODS, | Dickinson et al. 
20031) the C osmic Evolution Survey (COSMOS, | Sanders 
et al.|[2007 ), the Galaxy Mas s Assembly ultradeep Spec 



troscopuTSurvey (GMASS, Cimatti et al. 2008), the 



HST Cosmic Assembly Near-infrared Deep Extragalactic 



Legacy Survey (C ANDELS, |Grogin et at1|2011| |Koeke 



moer et al 
evolution o: 



20111), that have allowed us to study the 
galaxies from the earliest cosmic epochs 



However, a limitation of such surveys is the relatively 
small volumes prob e d, eve n at high redshifts: for ex- 



ample, Ilbert et al. (2006) find noticeable field-to-field 



variations in redshift distributions in the Canada- France- 
Hawaii Telescope Legacy Survey (CFHTLS 1 ) in fields of 
0.7- 0.9 deg 2 . 

Until lately, the combination of depth and area re- 
quired to map a large volume (~ 1 Gpc 3 ) of the high- 
redshift Universe at near-infrared wavelengths, where the 
redshifted emission from stars in distant galaxies peaks, 
has been prohibitively expensive in telescope time. Two 
recent developments have now made this regime accessi- 
ble. On the ground, the availability of wide-field near- 
infrared cameras has greatly improved the effectiveness 
of ground-based near- infrared surveys in the 1 — 2.5 /im 
wavelength range. In space, the exhaustion of the cryo- 
genic coolant of the Spitzer Space Telescope opened up 
an opportunity to pursue large near-IR surveys using 
the two shortest wavelength channels (IRAC1 [3. 6] and 
IRAC2 |"4.5 ] ) of the Infrared Array Camera (IRAC, [Fazio 
et al.|2004 ) in the post-cryogenic or "warm" mission that 
were much larger than was feasible during the cryogenic 
mission. The Spitzer Extragalactic Representative Vol- 
ume Survey (SERVS), a Spitzer "Exploration Science" 
program, stems from these two developments. 

SERVS is designed to open up a medium-depth, 
medium-area part of parameter space in the near-infrared 
(see Figure [l]), covering 18 deg 2 to ss 2 /xJy in the Spitzer 
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[3.6] and [4.5] bands. These observations required 1400hr 
of telescope time and covered five highly observed as- 
tronomical fields: ELAIS-N1 (hereafter EN1), ELAIS-S1 
(ESI), Lockman Hole (Lockman), Chandra Deep Field 
South (CDFS) and XMM-large-scale structure (XMM- 
LSS). The five SERVS fields are centered on or close 
to those of corresponding fields surveyed by the shal- 
lower Spi tzer Wide-area Infrar ed Extragalactic Survey 
(SWIRE; Lonsdale et al. 2003), and overlap with sev- 
eral other major surveys covering wavelengths from the 
X-ray to the radio. Of particular importance is near 
infrared data, as these allow accurate ph otometric red 
shifts to be obtained for high redshifts ( van Dokkum 
et~^[2006l IBrammer et al. 2008; jllbert et al |2TO[ [Car 
damone et al. 2010): SERVS overlaps exactly with the 
12deg z of the VISTA Deep Extragalactic Observations 
(VISTA VIDEO, Jarvis et al. 2012, in preparation) sur- 
vey (Z 7 Y, J, H and K s bands) in the south, and is cov- 
ered by the UKIRT Infrared Deep Sky Survey (UKIDSS 
DXS, 



Lawrence et al. 



2007) survey (J, K) in the north. 



SERVS also has good overlap with t he Herschel Multi 
tiered Extragalactic Survey (HerMES , |01iver et al.|2012 1 
in the far-infrared, which covers the SWIRE and other 
Spitzer survey fields, with deeper subfields within many 
of the SERVS fields. 

Sampling a volume of ~ 0.8 Gpc 3 from redshifts 1 to 
5, the survey is large enough to contain significant num- 
bers of rare objects, such as luminous quasars, ultralu- 
minous infrared galaxies (ULIRGs), radio galaxies and 
galaxy clusters, while still being deep en ough to find 
L* ga laxi es out to z ~ 5 (s ee for example Falder et al. 
|2011[ and Capak e t al.|20il"1 who find two galaxies in the 
Z = 5.3 cluster bright enough to be detected by SERVS 
at 4.5/um.) For comparison, the largest structures seen 
in the Mille nnium simulation at z ~ 1 are of the order 
of 100 Mpc ( |Springel et~alj|2005| ) , which subtends 3° at 
that redshift, so each SERVS field samples a wide range 
of environments. By combining the five different fields of 
SERVS, the survey effectively averages over large-scale 
structure, and presents a representative picture of the 
average properties of galaxies in the high redshift Uni- 
verse. 

Spitzer observations of the five SERVS fields are pre- 
sented in detail in Section[2j Image processing is detailed 
in Section[3j focusing on the mosaicking process and uni- 
formity ofcoverage. Section [4] presents the extracted 
SERVS catalogs, as well as an assessment of overall data 
quality, detection limits and expected number counts. 
Section [5] gives an overview of the ancillary data avail- 
able at aiffcrent wavelengths in the five fields. Prelim- 
inary science results and science goals are described in 
Section [6j A summary of the SERVS data at hand is 
providedin Section [7] 



2. SPITZER OBSERVATIONS 

2.1. Selection of fields 

SERVS consists of five fields located near the centers 
of corresponding SWIRE fields: EN1, ESI, Lockman, 
XMM-LSS and CDFS. The SWI RE fields are in region s 
with low infrared backgrounds (Lonsdale et al.||2003|), 
making them ideal for follow-up at far-infrared wave- 
lengths. The SERVS fields were selected to have good 
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Fig. 1. — Area versus depth for SERVS compared to other surveys at wavelengths of RJ 3.6/xm (left panel) and fti 4.5/^m (r ight panel). For 
consistency, the depth shown is the 5 a limiting flux for point sources, excluding confusion noise (<t pp as described in Section |4.2[ l, calculated 
from the Spitzer performance estimation tool (http://ssc.spitzer.caltech.edu/warmmission/propkit/pet/senspet /index. htmlj in each case. 
The surveys are (from left t o right): GOODS, the Spitzer rollow-up to the UANUiiLS HS'l survey (Uosmic Assembly Near-IR Deep 
Extragalactic Legacy Survey, |Grogin et al.|201l||Koekemoer et al.|2011[| , the Spitzer Extragalactic Deep Survey (SEDS, Program identifier 
- hereafter PID - 60022, 61040, 611)41, 61042, 61043, P.l. G. Fazio), the Spitzer IRAC/MUSYC Public Legacy in E-CDFS (SIMPLE) survey 
(Spitzer, PID 20 708), the Spitzer U ltra Deep Survey (SpUDS, PID 40021, P.I. J.S. Dunlop), S-COSMOS, the Spitzer Deep Wide-Field 
Survey (SDWFS, As hby et al.|2009|l , the Spitzer-KETDEX Exploratory Large Area (SHELA, PID 80100, P.I. C. Papovich) Survey, S WIRE, 
the SPT-Spitzer Deep Field (SSDF, PID 80096, P.I. S. Stanford) and the Wide-Field Infrared Explorer (WISE, [Wright et al.|2010| |. 



TABLE 1 
The geometry of the SERVS fields 



Field Name 


Field Center 


Field PA 


Field Area 


Vertices* 


of the area covered by both [3.6] & [4.5] 




RA, Dec (J2000) 


(deg) 


(deg 2 ) 




(deg) 




EN1 


16:10:00, +54:30 


350 


2.0 


(244.2,54.2) 


(243.1,55.4) (240.9,54.8) 


(241.7,53.6) 


ESI 


00:37:48, -44:00 





3 


(10.5,-44.9) 


(10.4,-42.9) (8.4,-43.0) 


(8.4,-45.1) 


Lockman 


10:49:12, +58:07 


328 


4.0 


(165.0,57.4) 


(161.7,59.8) (159.3,59.0) 


(162.7,56.4) 


CDFS 


03:32:19,-28:06 





4.5 


(54.4,-27.1) 


(51.8,-27.0) (51.7,-28.9) 


(54.4,-28.9) 


XMM-LSS 


02:20:00,-04:48 





4.5 


(37.2,-5.4) 


(37.0,-3.9) (33.9,-4.1) 


(34.3,-5.7) 


* Single-band catalogs extend beyond vertices. 



overlap with current and proposed surveys in other wave- 
bands within the SWIRE fields (see Section fl]), to cover 
both northern and southern hemispheres, ana to have a 
range in Right Ascension allowing both flexible follow-up 
with ground-based telescopes and good scheduling op- 
portunities for Spitzer. Field geometry and observation 
details are given in Table [j] & Table [2l The observed 
SERVS mosaics are shown in Figures ^2] through [6j to- 
gether with the coverage of significant overlapping sur- 
veys (see Section [5] for an exhaustive list of all ancillary 
data in and near the SERVS fields). 

A small fraction of the SERVS area was already cov- 
ered by other deep surveys with IRAC, such as the 
Spitzer IRAC /MUSYC Public Legacy in E- CDFS survey 



(SIM PLE, PID 20708; PI P. van Dokkum; |Damen et al . 
2009| in the CDFS field and the Spitzer Ultra Deep Sur- 
wyfSpUDS, PID 40021; PITS. Dunlop) in XMM-LSS. 
In order to minimize the total required telescope time, 
these particular areas were not observed. A selection of 



the IRAC [3.6] and [4.5] data from these surveys (both 
of which also use the 100 s frametime) are therefore sub- 
sequently added into the final SERVS mosaics to attain 
an app roxim ately uniform overall depth (see details in 
Section 3.3 1. 

In addition to this pre-existing Spitzer data, there are 
also two smaller deep fields located in the SERVS area: 
AORID 2 4402688 in Lockman (PID 64; P.I. Fazio) and 
the overlapping pointings of AORIDs 6005016 (PID 196, 
P.I. Dickinson) & 10092288 (PID 3407, P.I. Yan) in EN1. 
We deliberately reimaged them as part of SERVS so that 
the data taken during the post-cryogenic period could be 
compared with the data taken earlier in the mission, and 
their small size made tiling around them very inefficient. 

2.2. Design of observations 

2 An individual Spitzer observation sequence is an Astronomical 
Observation Request (AOR) 
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244.0 2-13.5 243.0 242.5 242.0 241.5 241.0 
Right ascension [deg. J2000) 



Fig. 2.— The [3.6] SERVS mosaic image of the EN1 field. Sur- 
veys of com parable size are sho wn here, such as HerMES level 5 
(in red, se e |01iver et al.||20I2] for details), CFHT H-band (dark 
blue), UKIDbb J &; tt.-bands (light blue), the Chandra X-ray sur- 
vey (green), SWIRE IRAC (dashed magenta) and SWIRE MIPS 
(dashed orange). Surveys encompassing the entire SERVS field 
such as HerMES Level 6 and the GMRT survey at 6f0 MHz are 
not shown here. More details about ancillary data coverage can be 
found in Section [5] 




166.0 165.0 164.0 163.0 162.0 161.0 160.0 159.0 

Right ascension (deg. J2000) 



Fig. 4. — The [3.6] SERVS mosaic image of the Lockman field. 
Only surveys of comparable sizes are shown here. Superposed are 
the HerMES Level 5 (in red) and Level 3 (magenta & orange), 
the Owen/Wilkes deep VLA (orange), the UKIDSS J, K coverage 
(cyan). The Chandra survey is displayed in green and the GMRT 
survey in blue. SWIRE IRAC & MIPS are shown as dashed dark 
magenta and dashed dark green. Surveys encompassing the entire 
field such as HerMES Level 5 are not shown here. More details 
about ancillary data coverage can be found in Section [5] 
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11.0 10.5 10.0 9.5 9.0 8.5 
Right ascension (deg. J2000) 



8.0 



7.5 



Fig. 3.— The [3.6] SERVS mosaic image of the ESI field. Sur- 
veys of comparable size are shown here, such as the HerMES- 
VIDEO field (in red), the ATCA/ATLAS radio survey (dark b lue), 
the VIDEO surve y (light blue), the deep ES1-XMM field of |Fer-| 
uglio ct al. (2008) (magenta), SWIRE IRAC (dashed orange) and 
bWlKJi MIPS (dashed green). Surveys encompassing the entire 
field such as HerMES Level 6 are not shown here. More details 
about ancillary data coverage can be found in Section [5] 
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54.S 54.0 53.5 53.0 52.5 52,0 51,5 

Righl ascension (deg, J2000) 

Fig. 5.— The [3.6] SERVS mosaic image of the CDFS field. 
Shown here are the HerMES Level 2 pointing (in red), the 
ECDFS/MUSYC survey (orange), the SIMPLE & GOODS sur- 
veys (magenta and green respectively). The VIDEO pointing is 
shown in cyan and the ATLAS radio survey in blue. SWIRE IRAC 
& MIPS are shown as dashed dark magenta and dashed dark green. 
Surveys encompassing the entire field such as HerMES Level 5 are 
not shown here. More details about ancillary data coverage can be 
found in Section [5] 
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37.5 37.0 36.5 36.0 35.5 35.0 34.5 

Right ascension (deg. J2000) 



34.0 



Fig. 6.— The [3.6] SERVS mosaic image of the XMM-LSS field. 
The two HerMES Level 4 & 3 pointings are shown in red and or- 
ange. The VIDEO pointing is featured in cyan and VVDS in ma- 
genta. The Eastern fields correspond to the SpuDS (blue), SEDS 
(green) and CANDELS (light grey). XMM-LSS is shown in (dark 
green) and extends beyond the SERVS CDFS field limits. SWIRE 
IRAC is shown as dashed green and SWIRE MIPS as dashed brown. 
Surveys encompassing the entire field such as HerMES Level 6 are 
not shown here. More details about ancillary data coverage can be 
found in Section [5] 



The design of the SERVS observations reflected several 
tradeoffs to ensure efficient use of the telescope, accurate 
filling of the fixed field geometries, and reasonably flex- 
ible scheduling. The SERVS depth was selected so that 
the confusion level just became significant; attempts to 
make it much deeper would require better ancillary data 
(e.g., GOODS) reaching in the confusion noise (the rate 
at which depth is achieved no longer decreases as the 
square root of exposure time. Within the constraints 
of the call for proposals, at this depth, SERVS is the 
largest area that can easily be surveyed and that has 
matching ancillary data. The depth of SERVS allows us 
to detect all massive (> lO n M0) galaxies out to z ~ 4 
(see section 4.3), essentially the entire range of redshift 
over which they are seen. SERVS can thus trace the 
evolution of these objects from their formation until the 
present epoch. As a consequence of those factors, the 
survey covers 18 deg 2 and reaches down to ~ 2 /iJy in 
the Spitzer [3.6] and [4.5] bands. 

Each field was observed in two distinct epochs, with the 
difference in time between the two epochs ranging from 
a few days to several months 3 . This allows us to reject 
asteroids, and also gives a better photometric accuracy 
by ensuring that most objects appear in very different 
places on the array in the two sets of observations. It 
arises from both a half-array offset in array coordinated 
between the two epochs and the fact that the time dif- 
ference between the execution of each epoch results in 

3 When the scheduling gap was in months, the AORs were re- 
designed at the time of the observations in order to maintain proper 
alignment of the tiles. 




Fig. 7. — A detail of the data coverage in the EN1 field at 3.6 fim. 
The image is 12" across. The cyan lines indicate the array edges 
of the individual input images and the greyscale underneath shows 
the mosaic coverage depth, ranging from 12 to 35 frames. The non- 
uniform depth of coverage s een above is reflected in the cumulative 
distribution function of Fig|12| 



a difference in the field rotation, and hence a different 
map grid for the two different epochs (Figure [7] shows 
the coverage of the two different epochs of observation in 
the EN1 field). 

Toward the end of the IRAC warm instrument charac- 
terization (IWIC), several tests were performed on vari- 
ations of the SERVS AORs to establish which observa- 
tion strategy was optimal. Three strategies were tested, 
all using the small cycling dither pattern, which allows 
for good coverage while ensuring that objects are shifted 
by a minimum of several arcseconds between observa- 
tions. The three strategies considered were: (1) two 
epochs of three dithered 200s frames, (2) two epochs 
of three dithered pairs of repeated 100s frames, and (3) 
two epochs of six 100s frames. In theory, strategy (1) 
is the most efficient and should result in a lower read- 
noise contribution. However, in practice, artifacts from 
bright stars were strong, reducing the effective area, and 
the radiation hit (hereafter rad hit) numbers were high 
(each array receives approximately 1.5 hits per second, 
each affecting on average two pixels, as detailed in the 
IRAC instrument handbook 4 ), resulting in a few rad hits 
leaking through into the final mosaic. There was also no 
measurable improvement in noise level compared with 
the other two options, which used 100s frames. Option 
(2) was almost as efficient as option (1), as only a frac- 
tion of a second was added to the overheads (the 200s 
frames have a longer readout time than the 100s frames), 
but image persistence effects were significant. Option (3) 
was therefore adopted, resulting in a very robust survey 
at the expense of only « 3% of extra observing time. 

The performance of the IRAC camera (both optically, 
in terms of PSF and array distortion, and in terms of 
sensitivity) was similar to cryogenic performance, hence 
the survey design was not modified because of array tem- 
perature changes. The sensitivity of the [3.6] band was 
affected at the 7% level by a change in the array bias be- 
tween the taking of the early and later SE RVS fields (see 
details of the calibration issues in Section 3.2 1, but this 



4 The IRAC instrument handbook can be found at 
http: // irsa. ipac . caltech.edu/data/SPITZER/docs/irac/ 
irac instrumenthandbook/ 
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TABLE 2 

Observing log and IRAC Instrument settings for the SERVS fields 



Hinlri p\ ti H f^r^r^fh 




IRAC 




Array T 


fa fil hias 

O • U \ LJ 1 Oil) 








(; m t~\ a i crn lei 


ID 


(K) 


CmVI 
V m v ) 


Cm VI 


ENl epoch 1 


2009-07-28 to 2009-08-01 


PCI 


61050 


31 


450 


450 


EN1 epoch 2 


2009-08-02 to 2009-08-05 


PCI 


61050 


31 


450 


450 


ESI epoch 1 


2009-08-06 to 2009-08-11 


PCI 


61051 


31 


450 


450 


ESI epoch 2 


2009-08-14 to 2009-08-18 


PC2 


61051 


« 29* 


450 


450 


Lockman epoch 1 


2009-12-10 to 2009-12-21 


PC10,PC11 


61053 


28.7 


500 


500 


Lockman epoch 2 


2009-12-25 to 2010-01-04 


PC11,PC12 


61053 


28.7 


500 


500 


CDFS epoch 1 


2010-11-01 to 2010-11-13 


PC33,PC34 


61052 


28.7 


500 


500 


CDFS epoch 2+ 


2009-10-13 to 2009-10-28 


PC6,PC7 


61052 


28.7 


500 


500 


XMM-LSS epoch 1 


2010-10-01 to 2010-10-18 


PC31,PC32 


60024 


28.7 


500 


500 


XMM-LSS epoch 2 


2011-02-19 to 2011-03-06 


PC42 


60024 


28.7 


500 


500 


* The array tcmpcrat 


ures were allowed to float in 


campaign PC2. 











Epoch 2 of CDFS was observed before epoch 1, which proved to be unschcdulablc in its originally planned slot. 



variation was not seen as significant enough to warrant 
a change in survey strategy. 

The mapping strategy used the small cycling dither 
pattern, which ensured full coverage with our map spac- 
ing of 280" . However each epoch of a SERVS field takes 
long enough to observe that the field rotation changes 
significantly between the start and end of a single epoch 
of observations. Therefore it needed to be robust against 
the 7° field rotation between AORs expected in an m 10 
day window in most SERVS fields. To allow for this, and 
to allow for accurate filling out of fixed field geometries, 
the SERVS AORs were kept relatively small (3x3 maps 
of 5' x 5' frames) and spaced closely enough to ensure 
overlap for the largest expected field rotation. The small 
AORs also had the advantage of being easier to schedule, 
allowing the placement of downlinks and the insertion of 
short non-SERVS observations. 

The total observing time for ENl, ESI, Lockman, 
CDFS and XMM-LSS were 153.4h, 231. 8h, 354.5h, 
319. lh 323. 2h, respectively. The mean integration time 
per pixel of the resulting SERVS mosaics is close to the 
design depth at « 1200s. There are, however, both re- 
gions of significantly deeper data where AORs and map 
dithers overlap, and shallower areas, particularly around 
the edges, or where one epoch is affected by scattered 
light from a field star. Uniformity of coverage of t he five 
SERVS fields is discussed in detail in Section [3~3l 



(SSC). These images have been dark-subtracted, flat- 
fielded, and have had astrometric and photometric cali- 
bration applied. A pipeline 7 originally used for process- 
ing SWIRE data was improved and applied to the SERVS 
data to further clean the frames of artifacts. Specifically, 
this pipeline fixed an artifact called "column pulldown" 8 
found near bright stars, and also corrected interframc 
bias offsets by setting the background equal to that of 
a COBE-based model of the zodiacal background (the 
dominant background at these wavelengths). Due to the 
inability to use the IRAC shutter 2 , all IRAC data suf- 
fer from a variable instrument bias level known as the 
"first-frame effect" which varies depending on the recent 
detector history. Thus no measurement of the true in- 
frared background, nor of any spatial structure within 
the background larger than the array size of 5', is possi- 
ble. 




3. DATA PROCESSING 

SERVS data are available from the Spitzer Heritage 
Archive 5 . SERVS mosaics (image, coverage, uncertainty 
and mask mosaics) and catalogs, including ancillary data 
at other wavelengths taken as part of SERVS will be 
made available to the community during the summer of 
2012, ultimately through the Infrared Science Archive 
(IRSA) . Catalogs containing the full data set of ancillary 
data will be described in detail in Vaccari et al. 2012., 
in preparation 

3.1. Image post-processing and mosaics 

Data processing begins with the Basic Calibrated Data 
(BCD) image, produced by the Spitzer Science Center 6 

5 http : / /sha. ipac . caltech.edu/applications/Spitzer/SHA/ 

6 http : 77ssc . spitzer . caltech.edu/! 



Fig. 8. — Left: An example of a SERVS reduced single 100s 
exposure frame at [3.6]. The frame size is ~ 5' X 5'. Right: the 
final co-added [3.6] image. 

The data were co-added (see Figure [8]) using the 
Mopex 9 package available from the SSC (parameters 
used are listed in AppendixjB]). All the data from a single 
field were co- added onto a single frame; the two different 
wavelengths are reprojected to the same astrometric pro- 
jection so that their pixels align one-to-one. The data are 



7 http : //irsa. ipac . caltech . edu/data/SPITZER/SWIRE/ 

8 The column pulldown effect, which manifests in the slow read 
direction (columns) of the detectors at 3.6 and 4.5/xm, is a depres- 
sion in the zero-level of the column. 

9 Mopex and its associated documentation can be ob- 
tained at http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 
dataanalysistools/tools/mopex/ 
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reprojected with a linear interpolation onto a pixel scale 
of 0.6", providing marginal sampling at 3.6 /im. The mul- 
tiple dithers allow at least some recovery from the severe 
undersampling of the IRAC camera at these wavelengths. 

The depth reached by the SERVS observations can eas- 
ily be put in perspective when comparing SWIRE and 
SERVS cutouts of a similar region of sky (left and right 
respectively in Figurej9|. With a depth of 3.7 //Jy at [3.6] 
and 7.4 ^Jy at [4.5], SWIRE is limited toz« 1.5 for L* 
galaxies, whereas SERVS can detect these galaxies up to 
z ~ 4 (see also Figure 17 1. 




Fig. 9.— A cutout of the EN1 field at [3.6] (top panel) & [4.5] 
(bottom panel), as imaged by SWIRE (left) k, SERVS (right). The 
difference in depth between the two surveys can clearly be seen 
here. 

Large mosaics containing all the SERVS epochs of ob- 
servation in each field were created from the co-added 
images, along with th e uncertainty, coverage and mask 
images (see Figure 10 1. The resulting mosaics cover ar- 
eas of ~ 2, 3, 4, 4.5 and 4.5 deg 2 for EN1, ESI, Lockman, 
CDFS and XMM-LSS, respectively (see Table [I] for de- 
tails about the geometry of the five fields). 

3.2. Calibration issues 

The initial observations (all of EN1 and the first half 
of ESI) were made before the IRAC detectors were sta- 
bilized at their current operating temperature of 28. 7K, 
which is now the norm for the Spitzer post-cryogenic (i.e. 
"warm") mission (see Table [2]). Instead, these data were 
taken at a controlled temperature of 31K. The second 
half of ESI was taken during a time when the detec- 
tor was cooling from 3 IK to the final temperature of 
28. 7K, and was not under active temperature control. 
In addition, during this time period the detector biases 
for both arrays were adjusted. These changes in tem- 
perature and detector operating parameters resulted in 
changes to both the instrument calibration and its noise 
properties. These were only measurable in the [3.6] band, 
resulting in an 7% increase in the noise in EN1 and 
ESI at [3.6] compared with the remainder of the survey. 

IRAC is calibrated based on dedicated calibration ob- 
servations collected during science operations. These 



data are typically collected on timescales of years. The 
Spitzer cryogenic mission, as well as the nominal warm 
mission, are extremely well-calibrated. However, during 
this transition period only a small amount of calibration 
data could be taken, as conditions were constantly chang- 
ing. The Spitzer Science Center provided an initial cali- 
bration, which was reliable to a few percent. Photometry 
from SERVS was compared with that of SWIRE on an 
object-by-objcct basis, and small multiplicative offsets 
at the few percent level were found. These calibration 
errors were fixed at the catalog level (see Section [4]) by 
applying multiplicative factors derived from comparison 
between sources detected both in SWIRE and SERVS. 
Calibration is therefore the same as that of the SWIRE 
data. 

3.3. Uniformity of coverage 

Certain areas of the original SERVS footprint were al- 
ready covered by several previous Spitzer surveys, such as 
SpUDS for XMM-LSS and SIMPLE/GOODS for CDFS 
(see Section 2.1 for details). These regions were deliber- 
ately avoided during the SERVS observation campaign to 
save observation time. Since the optical properties of the 
camera did not change between the cryogenic and warm 
missions and since all the Spitzer data was combined us- 
ing the same Mopex pipeline, previous survey frames 
were subsequently merged with the SERVS data at co- 
addition and patched onto the final mosaics. Given the 
higher depth of the previously existing smaller Spitzer 
surveys, it was always possible to reasonably match the 
depths of the SERVS observations in these areas by se- 
lecting the right number of single frames to co-add. In 
addition, since the archival data are only a minimal frac- 
tion of the survey (~ 12 % of XMM-LSS and - 3 % of 
CDFS, none in the other fields), these minor differences 
in coverage between archival and SERVS data are not an 
issue. 

The selection of archival data was focused on filling 
the missing areas as best as possible and thus not refined 
to be necessarily chosen from different dates of observa- 
tions, hence some contamination by transient sources on 
timescales of hours is possible. Some artifacts such as 
muxbleed 10 were present in the cryogenic data and not 
in SERVS, however after processing, all artifacts were 
removed and did not impact the mosaicking. A grid 
extending over the missing data was set to match the 
SERVS frame centers (spacing and orientation) as closely 
as possible. The closest SpUDS (or SIMPLE/GOODS) 
frames to these grid centers were then automatically se- 
lected and a subsample of those was chosen to best cover 
the area and match the SERVS depth (the left panel of 
Figure 11 shows the XMM-LSS field frame selection pro- 
cess as an example). The resulting coverage map for the 
XMM-LSS field is shown in the right panel of Figure [Tl] 
as an example. 

Coverage is thus not completely uniform throughout 
the fields, though averagin g ~ 1400 s of exposure time 



over all five fields (see Figure 12 ). By design some overlap 
between AORs was allowed and some areas in a field may 
have a higher coverage (e.g. at the AOR intersections), 

10 The muxbleed effect appears as a series of bright pixels along 
the fast read direction (horizontal in array coordinates), and may 
extend the entire width of the array (5'). 
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Fig. 10. — SERVS final set of mosaics for the EN1 field at 3.6/mti (top left: image mosaic, top right: coverage mosaic, bottom left: 
uncertainty mosaic, bottom right: mask mosaic). 



for example totalling up to ~ 4700 s in Lockman. Al- 
though great care was taken to minimize variations in the 
coverage throughout a field, the exposure time can rise up 
to ~ 5000 s (e.g. in XMM-LSS) when the added ancillary 
data intersects both SERVS epochs (mostly around the 
edges of the SERVS / ancillary data, as is obvious in the 
right panel of Figure 11 ). However these areas of higher 



coverage represent a very small portion of the fields and 
do not have an impact on the overall uniformity, given 
the already intrinsic nonuniform nature of the SERVS 
coverage due to the original tiling design. The higher 
depths of the observations at the frames overlaps or due 
to archival data overlapping the SERVS frames do not 
result in any impact to the survey reliability or com- 
pleteness when a minimum flux selection based on the 
SERVS data is used. 

3.4. Image masks 

Bright source artifacts (straylight, column pulldown, 
latents) are a significant problem with Spitzer data. Most 
of these artifacts are discussed in more detail in the IRAC 
instrument handbook 4 , available from the SSC, as well 



TABLE 3 

Lookup table for the mask image bit values 



Bit number Bit value 



Flag 



BIT00 = 1 Overall data quality 

BIT01 = 2 Set if pixel contains rad hit 

BIT02 = 4 Set if optical ghost present 

BIT03 = 8 Set if stray light present 

BIT04 = 16 Set if saturation donut 

BIT05 = 32 Set if pixel contains latent image 

BIT06 = 64 Set if pixel is saturated 

BIT07 = 128 Set if column pulldown present 

BIT08 = 512 Set if bright star is present 



Notes: Flags arc described in more detail in the text in Section 3.4 



as in Surace et al. (2004) and Surace et al. (2005 1 7 . In 
order to flag tor these artifacts, mask images are created 
for each of the five fields. Flags bit values attributed to 
each significant artifacts are listed in Table [3j as well as 
in the headers of the mask FITS files. 

In addition to these known Spitzer image artifacts, 
saturated stars are common across the deep and wide 
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Fig. 11. — Left: XMM-LSS single frame centers, showing epoch 1 in cyan and epoch 2 in blue crosses (see Table[2]l. A grid extrapolating 
the XMM-LSS epoch2 geometry was overlaid onto the SpUDS field frames (in green open circles) and used to select the relevant SpUDS 
frames to co-add (magenta solid points) to complete the XMM-LSS field. Right: Coverage map in equatorial coordinates of the XMM-LSS 
field at [3.6]. Note the addition of the SpUDS data between 35 < a < 34 & -5.7 < 6 < -4.7 to patch the XMM-LSS field. The rotation 
in tiling induced by the choice of two different epochs of observation (of six 100 s frames each, as discussed in Section |2,2| is also clearly 
visible here. Coverage is not completely uniform through out the field and data can be deeper, especially at the overlap of the SERVS 
(warm mission) & SpUDS (cryogenic) surveys (see Section [33] for details). 




EN I 

ESI 

LCKMN 

CDFS 

XMM 




1000 2000 3000 

Exposure (seconds) 



A 000 



Fig. 12. — Cumulative distribution function of the coverage maps 
for all five fields, in terms of exposure time (in seconds). The curve 
represents the fraction of pixels with that integration time or lesser. 

SERVS fields and thus need to be nagged appropriately. 
Indeed very bright stars do not have reliably extracted 
fluxes in SERVS, and may occasionally be split into mul- 
tiple fainter sources, or are saturated, so no flux can be 
measured accurately; erroneous detections in the wings 
of the PSF can also cause artifacts. As a result, a safe ra- 
dius has to be set and flagged around those bright stars. 
Luckily, any object triggering bright star artifacts in the 
SERVS data is e asily detected by the T wo Micron All Sky 
Survey (2MASS |Skrutskie et al.|2006h , which has reliable 
fluxes even for very bright objects. The positions and K- 
band magnitudes of bright stars are downloaded from 
the 2MASS Point Source Catalog catalog in the Vizier 
database 11 . Magnitudes are thus converted into radii ac- 



FlG. 13. — Bright star mask cutout of a portion of the Lockman 
field. Black circles show the masked regions around the cross- 
identified 2MASS bright stars centers (red crosses) in the field. 
Radii are proportional to 2MASS K-band brightness and shown in 
Table H 



cording to Table [4 The radii are taken from SWIRE 
(Surace et al. 2004) and increased to take into account 
the deeper exposures of SERVS. To reduce memory re- 
quirements on such large images, subsets of the mosaics 
were cut around each bright star, flagged and then reem- 
bedded onto the final mask mosaic. The resulting bright 
star flag masks a circular region around each bright star 
in the fields. A part of the Lockman field bright star 
mask is shown in Figure |13| 



4. CATALOGS 
4.1. Catalogs extraction and calibration 



The Vizier database is available at http://vizier.u-strasbg.fr/ 



Catalogs we re made using SExtractor (Bertin & 
Arnouts 1996). Two sets of catalogs were produced tor 
each field (the parameters used in the source extraction 
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TABLE 4 

Bright Object flagging for the mask images 



2MASS-K mag. range 


radius (") 


> 12.0 





10.0 - 12.0 


15 


8.0 - 10.0 


20 


7.5 - 8.0 


30 


6.5 - 7.5 


45 


5.0 - 6.5 


60 


< 5.0 


120 



Notes'. The radii arc taken from SWIRE (S uracc et al.||2004| l and in- 
creased to take into account the deeper exposures of SKHVy. 



TABLE 5 
Aperture sizes and corrections 



Aperture 
number 


Aperture radii 
(arcsec) 


[3.6] correction 
(arcsec) 


[4.5] correction 
(arcsec) 


apl 


1.4 


0.585 


0.569 


ap2 


1.9 


0.736 


0.716 


ap3 


2.9 


0.87 


0.87 


ap4 


4.1 


0.92 


0.905 


ap5 


5.8 


0.96 


0.95 



Notes: Aper ture sizes and corrections were derived for SWIRE |Surace| 
ct al. (2004). More details can be found therein. 



are given in AppendixlB I . The first set is based on extrac- 
tions from the [3.6] and" 4.5] bands, separately. A second 
set is then created when the two catalogs are merged, and 
only detections common to both are retained. This re- 
sults in a high-reliability catalog that is used to combine 
with other data sets. For rare objects searches, we rec- 
ommend using the high-reliability catalogs. The single- 
band catalog may be used for bandmerging with other 
data (e.g. UKIDSS/DXS data), or for statistical studies. 

SExtractor aperture fluxes are computed within radii 
of 1.4", 1.9", 2.9", 4.1" and 5.8" and corrected using 
the ap erture correction factors derived for SWIRE 
DR2/3( |Surace et al.||2004[ ) and reported in Table [| The 
IRAC instrument has a calibration tied to standard stars 
as measured in a fiducial 12" radius aperture. This aper- 
ture is non-ideal for faint-source extraction, so the fluxes 
are measured in smaller apertures and suitable aperture 
corrections are applied. 

A fraction of the SERVS data was taken prior to the 
final temperature stabilization of IRAC, and prior to the 
selection of the final array biases (affecting the EN1 and 
ESI fields, see Table [2]). During this "floating tempera- 
ture" period, the IRA(J calibration drifted, with changes 
to both the overa ll ga in and the detector linearization, as 
discussed Section 3.2 As a result, the images have overall 
calibration errors at the level of a few percent. In addi- 
tion, the other fields show small, but noticeable calibra- 
tion differences compared to the cryogenic [3.6] SWIRE 
data. These calibration errors were fixed at the catalog 
level by applying multiplicative factors derived from com- 
parison between sources detected both in SWIRE and 
SERVS (/ [3 .6]swire//[3. 6 ]servs = 1-07 for EN1 & ESI 
and 1.02 for Lockman, CDFS and XMM-LSS; at [4.5] 
the correction factors are very close to unity, and no cor- 
rections were applied) . SERVS calibration is thus tied to 
that of the SWIRE data. 

4.2. Number counts and completeness 

Simple number counts were derived from the extracted 
catalogs. The SERVS source counts for the XMM- 
LSS fie ld is provided as an example in the left part of 
Figure 14 (the remaining four fields are shown in Ap- 
pendix XF. Several features are visible in this plot. Fun- 
damentally, the observed counts present as a power-law. 
This is shallow at the bright end, primarily due to the 
presence of bright stars. There is a known break in the 



power-law index near 100 /LtJy (Glazebrook et al. 1994) 



At the faintest levels, the turnover is a result of the in- 
creasing incompleteness of the survey. It is clear that the 



SERVS completeness level catastrophically drops near 2- 
3 //Jy. As a consistency check, a sample of bright stars 
(SERVS fluxes within 0.3 < / < 2mJy and stellar- 
ity index > 0.95) was selected. Cross-identifying with 
the 2MASS catalog, K-band measurements were used 
to compare the deviation of SERVS colors with respect 
to the 2MASS ones. A color-color plot of i^MASS — 
M 3 . 6 , im (Vega) versus M 3 . 6/ _ tm (Vega) - M 4 . 5Mm (Vega) for 
the XMM-LSS field, shown in the right part of Figure [Til 
confirms that SERVS is consistent with 2MASS. 

SERVS detects ~ 100,000 sources per square degree 
and, with ~ 40 beams per source, approaches the clas- 
sical definition of where source confusion becomes im- 
portant. The effective depth of SERVS is thus affected 
by confusion noise, which makes the definition of survey 
depth dependent on the experiment one wishes to carry 
out. For a detection experiment on point-source objects 
detected in another waveband with a well-determined 
position (positional uncertainty much lower than the 
SERVS beam of ps 2"), the appropriate number, denoted 
CTpp, is determined from the pixel-to-pixel variance in a 
source-free region of a single basic calibrated data (as- 
suming it scales as the square root of the coverage and 
that extraction is carried out by source fitting; i.e. as- 
suming 7.0 and 7.2 noise pixels in [3.6] and [4.5], respec- 
tively). Another measure of the noise, which includes 
some contribution from confusion noise, is obtained from 
"empty-aperture" measurements, where the standard de- 
viation of fluxes in object-free apertures in the final mo- 
saic is measured directly. For our fields, this measure- 
ment was made in 3.8" diameter apertures (SWIRE aper- 
ture 2), and is denoted cr ap . SWIRE aperture 2 is rec- 
ommended by the SWIRE team as the most stable aper- 
ture for photometry since most faint IRAC sources are 
slightly resolved at the ~ 1 — 2" level. Survey depths for 
the various measurements are summarized in Table [6l 

Finally, for survey work, the completeness limits at 
50 and 80% (S C 5o and S C 8o)j give a good indication of 
the depths that are usable for global survey pro perties 
such as sour ce c ounts. Following t he techniques of |Chary] 
et^p004] ) and |Lacy et al^|2005l ), we simulated -10,000 
model galaxies, distributed them in the reduced mosaics, 
and extracted them using the same pipeline that cre- 
ated the catalogs. The recovery rate of the model galax- 
ies was used as our completeness indicator and suggests 
that « 50% completeness is reached at a flux density of 
~ 2 — 3 /iJy in the single band catalogs at both [3.6] and 
[4.5], due to a combination of signal-to-noise ratio and 
source confusion (see Figure 15 1. More details of the 
completeness and reliability oftlie SERVS catalogs will 
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1358 [3.6] sources 
0.3 mjy < F < 2 mJv 
&nd stellarily > 0.95 
with 2Mass K-band . 




Fig. 14. — Top: SERVS number counts versus flux at [3.6] 
(black histogram) and [4.5] (red histogram) for the XMM-LSS field. 
Grey dashed lines show the selection of sources used in the right 
plot (with fluxes as 0.3 < / < 2mjy). Bottom: Color-color 
plot showing _fsT 2 MASS - M3.6 M m(Vega) versus M3.6 Mm (Vega) - 
M A . (Vega) for sources within the flux range defined by the 
grey dashed lines above, plus a cut in stellarity index > 0.95 and 
the existence of a 2MASS K-band measurement as an additional 
constraint. Red dashed lines help pinpoint the location of the (0,0) 
point in this diagram. Similar plots for the EN1, ESI, Lockman & 
CFDS fields are shown in Appendix [A] 



be presented by Vaccari et al. (2012, in preparation). 

4.3. Expected detection limits & redshift distribution 

The SERVS project uses semi-analytic models exten- 
sively, both to make testable predictions of the properties 
of SERVS galaxies, and to inform our follow-up strate- 
gies in wavebands other than the near-infrared (e.g. Fig- 



18 1. Two different semi-analytic models are currently 



being used. 

The first set is based on the Guo et al. (20111 ver 



TABLE 6 

Approximate survey depth and completeness in the SERVS 

FIELDS 



Measurement 


[3.6] value 


[4.5] value 




(MJy) 


(^Jy) 


5<Tpp 


1.3 


1.5 


5 (Tap 


1.9 


2.2 




4.0, 3.0 


3.5, 3.5 


7, 5 


5, 5 



t The two completeness levels shown correspond to the the two-band 
detection catalog and the single-band catalog. 

Notes: EN1 and ESI are about 7% noisier than the other fields in 
[3.6] due to different detector settings used early in the warm mission 
(Table 
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Fig. 



15. — Completeness plots for Lockman Hole at [3.6] shown 
in green, [4.5] shown in blue, and the reliable dual-band detection 
cat alogs shown i n blac k. Using the techniques of ]Chary et al.| | [2004[l 
andlLacy et al.l (120051 , ~10,000 model galaxies were simulated and 
placed m the mosaics. The sources were extracted using the same 
pipeline as the catalogs and the completeness was estimated as the 
recovery rate of the simulated model galaxies. 



cones 12 were created and fully described in Henriques et 
al. (2012). The lightcones contain a wide range of pho- 
tometric bands that cover the UV to near-infrared region 
of the spectra and allow an exact match to the observed 
selection criteria. They a lso include a choice betw een 
flu xes computed usi ng the Bruzual & Chariot ( 2003 1 or 
the Maraston ( 2005 ) stellar populations. The latest have 
been shown to reconcile the predicted X-band rest-frame 
luminosity function w ith observations at high redshift 



(Henriques et al. 2011 Henriques et al. 2012), for which 
1RAC data - such as those obtained with SERVS - have 
been essential. The second set, from van Kampen et 
al. (2012a, in preparation), includes both the effects of 
halo-halo and ga laxy-galaxy mergers, and uses GRASIL 
( |Silva et al.|1998 ) to predict spectral energy distributions 
(SEDs) from the optical to submillimeter. 

The expected redshift distribut ion is derived from the 
simulations and shown in Figure |16| A lower flux limit 
of 2 piJy (AB magnitude of m [3 6] < 23.1) for the SERVS 



sion of the Munich semi-analytic model for which light 



12 These lightcones are pu blicly available at |http://www.mpa-| 
garching.mpg.de/millennium 
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— Henriques2012 - M05 (l tJAm) <23.1) 
-- Henriques2012 - BC03 (l (3J( „ ) <23.1) 



Fig. 16. — SERVS expected normalized galaxy counts at the de- 
tection limit of 2 (ijy (or mra gi < 23.1). The two histograms cor- 
respond to flux lim its applied to the [3 . 6] ban d either using the 
|Maraston] | |2005[ l or JBruzual fc Charlot| ( |2003| ) stellar population 
models (see Section |4.3|. 
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survey was used here. The expected redshift distribution 
of the SERVS galaxies peaks around z ~ 1 and extends 
to redshifts of z ~ 3, with a small fraction of objects 
detected up to z ~ 4 (for basic comparison purposes, the 
SWIR E photometr ic redshift distribution can be found 
in [Rowan- Robinson et al. 2008; at the current time, the 
SERVS photometric redshift distribution has not been 
derived but will be presented by Pforret al. 2012b, in 
preparation, as discussed in Section 5.7). Galaxies de- 



tected in SERVS will therefore span the epochs where 
galax ies gain the yast ma jori ty of their stellar mass. In- 



deed, prownetaL (J2007J) and Cool et al. (2008) estimate 
that L" galaxies roug hly double in mass betwe en z = 
and 2Ril. In addition, | van Dokkum et ah ( 2010 ) recently 
showed that about half the mass ot any given large galaxy 
is added between z — and z = 2, by comparing galaxy 
samples at constant number densities. SERVS will be 
able to extend such studies out to higher redshifts with 
good statistics. 

In order to roughly estimate the luminosity of the 
faintest galaxy SERVS will be able to detect in the IRAC 
[3.6] band as a function of redshift, three model SEDs 
were used: a starburst similar to M82, a 250 Myr stellar 
population and a 5 G yr stellar population from |Maras-| 



ton 



( 2005 1 . Figure [l7| shows the luminosity of the faintest 
detectable source from 0.1 fim to 1mm as a function of 
redshift, if this source was represented by one of the three 
SED models we used. SERVS might be able to detect L* 
galaxies up to z ~ 4 and 0.1L* out to z ~ 1. The lu- 
minosities of the Maraston SEDs can easily be converted 
to stellar masses by multiplying the luminosity of the 5 
Gyr and 250 Myr stellar population models by 1.15865 
and 0.00133 respectively. SERVS will ensure the deriva- 
tion of robust stellar masses because it includes the rest- 
frame near-infrared out to high redshifts. This coverage 
is essential to break the degenera cy between star forma - 
tion history and dust reddening (Maraston et al. 2006). 
Photometric redshif ts fo r SERVS galaxies are also being 
derived (see section 5.7). 



Fig. 17. — Top panel: SERVS faintest galaxy to be detected, in 
units of solar luminosities (integrated luminosity from the optical 
to the far-IR: 0.1 /zm— 1mm) at a given redshift. Three Spectral 
Energy Distribution (SED) models are considered here and shown 
in the bottom panel. The three SEDs consist of an M82-like star- 
burst galaxy (in blue), a 250 Myr stellar population (magenta) and 
a 5 Gyr stellar population (red), scaled to the same 3.6/^m rest- 
frame. SERVS can detect L* galaxies out to z ~ 4, and 0.1L* out 
to z ~ 1. At redshift 5, a starburst galaxy of luminosity ~ 10 12 Lq 
could potentially be detected. 



5. ANCILLARY DATA 

Overlaps with existing and surveys in progress are 
listed in Table [7] Here we provide a more detailed de- 
scription of some of the most significant (in terms of their 
overlaps with SERVS) of these surveys. 

5.1. Optical surveys 

A number of optical surveys overlap one or more of 
the SERVS fields. Among the most signific ant are the 
ESQ /Spitze r Imaging Survey (ESIS) in ESI ( |Berta et al 



20061 120081) which reaches depths of 25, 25, 24.5, 23.2 
(Vega) in B, V,RSzI, respectively, the Canada-France- 
Hawaii Telescope Legacy Survey (CFHTLS; both deep 
and wide pointings in u*,g,r,i & z in the XMM-LSS 
field reaching depths of 28.7, 28.9, 28.5, 28.4, 27.0 and 
26.4, 26.6, 25.9, 25.5, 24.8 AB magnitudes, respectively), 
and ancillary SWIRE data in Lockman, EN1 and CDFS 
in a variety of de pths and filters, but typically reaching 
at least r — 24.5 (jGonzalez-Solares et aL|l2011). 



The Spitzer Adaptat ion ot the Red Sequence Cluster 
Surve y (SpARCS 13 , |Muzzin et aT " 



2009 Wilson et al. 



20091 DeGroot et al., 2012, in prep, Muzzin et al., 2012a 



in prep) has imaged the entire SWIRE area (excluding 
the XMM-LSS field which is covered by the CFHTLS) 
in the z' filter using MegaCam on the Canada France 
Hawaii Telescope (CHFT) or the Mosaic camera on the 
Blanco Telescope at Cerro Tololo Inter-American Ob- 
servatory (CTIO). The MegaCam observations reach a 

13 http:/ /www. faculty.ucr.edu/~gillianw/SpARCS 
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TABLE 7 

Surveys at other wavelengths covering >10% of a SERVS field (all-sky surveys excepted). 



Field X-ray Optical Ncar-IR Mid-IR Far-IR/ Radio 

Name data Data Data Data submm Data Data 



EN1 


N04 1 


SWIRE/INT 2 
SDSS 3 
SpARCS 4 


DXS 5 


SWIRE 6 


HerMES 7 L5 
S2CLS 8 


G08a 9 
GrlO 10 


ESI 


F08 11 


ESIS 12 
VOICE 13 , SpARCS 4 


VIDEO 


SWIRE" 


HerMES L6 


ATLAS 


Lockman 


W09 16 


SWIRE/KPNO+INT 2 
SDSS 3 
SpARCS 4 


DXS 5 


SWIRE 6 


HerMES 7 L3,L5 
S2CLS 8 


OM08 17 
G08b 18 ,G10 19 
I09 20 


CDFS 


CDFS 21 


SWIRE/CTIO 22 


VIDEO 14 


SWIRE 6 


HerMES 7 L2,L5 


ATLAS 15 






GaBoDS 23 
VVDS 24 




SIMPLE 25 


LABOCA/LESS 26 


M08 27 






VOICE 13 , SpARCS 4 










XMM-LSS 


XMM-LSS 28 


CFHTLS 30 


VVDS 31 


SWIRE 6 


HerMES 7 L5 


VVDS-VLA 32 




SXDS 29 


VVDS 24 
SXDS 34 


UDS 5 
VIDEO 14 
DXS 5 


SpUDS 25 


S2CLS 8 


S06 33 



Notes: [1] Chandra proposal 6900602 CP.I. N a ndral: [2] |Gonza,le,z-Solares et al.J (2011}: [31 |Abazaiian_ et al7| l|2009): [41 DeGro ot et 
al. 2012, in preparation [ 5] |Lawrence et al.| l|2007l; [61 |Lonsdale et al.| j2003| l, THurace et al.| j^UU5| >~- |7| | oliver"et all | 2012l; [8] 
http://www.jach hawaii.edu/JCffIVsurvcys/CosmoloKy.html [9| |Garn et al.|l |20u8a| l; 1 10] Grant ct al.| ]2010[ ~|Tl] |Feruglio et al.|]2008[ ; [12 |[Berta 
ct al. (2006 2008); 1 13 1 P.l. G. Covonc He M. Vaccari, http://pcoplc.na. inm.it/^covonc/voicc/ voicc.html |14| Jarvis ct al. (2012), in prep; [15] Norris 
et al.| l [27?U^ , |Middelberg et al.|||2008|l; [16l| Wilkes et al.| l |2009t ; |17| |Owen Hz Morrison| l |2008t ; ]18] |Garn et al.|l|2008b( ; [19] |Garn et al.fooioj; [20||lbar 
et^^ J2TlU9 t ;^^TT^eTrme^Tt^ujl2' Hu ? 5} ; [22]" http:/ / www.ast ro.caltech.edu/~bsiana/cdis_opt [23J G arching- Bochum Deep yurvey, pilc!e^randteT^T 
j2Uu6[ ; [24[ |Le Fevre et al.|]2005[ ; [25] |Damen et al.j ( |2009[ > , irsa.ipac.caltcch.edu/data/yPlTZER/docs/spitzcrmission/obscrvingprograms/lcgacy/ ; 
|26 yurvey ot the GUi y with the Large Apex Bolometer Camera (LAB OCA),|Weifi et al. 102009); [271 IMiller et al.|j2008 |; [28] |Pierre et al.M2007t ; 
[29] lUeda et al.l (|2008 [>; [3 0] www.cfht.hawaii.edu/ Science/CFHTLS; [31] |lovino et al.| ]20o5[ ; |Temporin et al.| l |2u08[ ; [32] |ijondl et al.| l |2u07j ;"[33] 
|Simpson et al.|(2006f ; [34] |Furusawa et al.| l |2008[ ; 



mean depth of 24.2 AB magnitudes, and the Mosaic cam- 
era observations reach a mean depth of 24.0 AB magni- 
tudes. The SpARCS collaboration has spectroscopically 



confirmed ~ 15 clusters at z > 1 (Muzzin et al. 
Wilson et al."1[2009l |Demarco et ah] jj(JlU[ |Mi~ 



2009, 
et al. 



uzzm 

2011[ Muzzin et al.', 2012b, in prep, Wilson et al., 2012a, 
in prep), and is in the process of carrying out detailed 
multipassband and spectroscopic follow-up studies (Ret- 
tura et al., 2012, in prep, Lidman et al., 2012, in prep, 
Noble et al., 2012, in prep, Ellingson et al., 2012, in prep, 
Wilson et al., 2012b, in prep). 

Started in October 2011, the VST/VOICE survey (P.I. 
G. Covone & M. Vaccari) is surveying the CDFS and ESI 
fields in u, g, r, i aiming at reaching AB^ 26 at 5 a. 

Optical spectroscopy has thus far been confined to 
small regions of SERVS - such as the re cently completed 
PRIMUS (PRIsm MUlti-object Survey , |Coil et al.|2011 | 
survey, covering parts of the ESI, CDFS & XMM-LSS 
fields - or to specific types of objects. The largest spec- 
troscopic survey is the VVDS, which has ps 10 000 spec- 
troscopic redshifts for field galaxies with 1 7.5 < Jab < 24 
in th eir XMM-LSS and CDFS subfields ( |Le Fevre et~aT 
2005 1 . In addition, spectra of AGN (Active Galactic Nu- 
clei) and quasars, now totaling several hundred objects 
selected by vario us techniques, hav e been obtained in the 



Trichas et al. (2010) 



SERVS fields by Lacy et al.| (|2007|) 
and Lacy et al. (2012, in prep 



5.2. Ground-based near-infrared surveys 

One of several major surveys to be carried out by 
the Visible and Infrared Survey Telescope for Astron- 



omy (VISTA) is the five filter near-infrared VISTA Deep 
Extragalactic Observations (VIDEO) survey (Jarvis et 
al. 2012, in preparation), which will cover 12 deg 2 to AB 
magnitudes of 25.7, 24.6, 24.5, 24.0 and 23.5 in Z, Y, J, H 
and K s filters. The ESI, XMM-LSS and CDFS SERVS 
fields are designed to exactly overlap their corresponding 
VIDEO fields. The combination of VIDEO and SERVS 
will be a particularly potent tool for the study of galaxy 
evolution at high redshifts. 

The Deep eXtragalactic Survey (DXS) is p art of the 
UKIRT Inf rared Deep Sky Survey (UKIDSS , |Lawrence 
et al.|[2007 ), and will cover the Lockman and EJN1 fields 
to 23.1 and 22.5 (AB) in J and K respectively. As of 
October 2011, Data Release 7 is available to the entire 
astronomical community. It covers parts of the EN1, 
Lockman and XMM-LSS fields. Data Release 9 is avail- 
able to the community served by the European Southern 
Observatory, and includes additional data. 

5.3. Mid- and far-infrared, and submillimeter surveys 

The SERVS fields were designed to be contained within 
the SWIRE fields. This has mostly been achieved, al- 
though constraints from other surveys mean that a small 
fraction of SERVS lies outside of the SWIRE coverage. 

The HerMES 14 survey is a Herschel Key Project to 
survey most of the extragalactic Spitzer fields, including 
the SWIRE fields. HerMES has six levels, correspond- 
ing to increasing depths, level 6 being the shallowest. 
Smith et al. (2012 1 give measured flux densities at which 



http:/ /hermes. Sussex. ac.uk/ 
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50% of injected sources result in good detections at the 
SPIRE wavelengths of (250, 350, 500) /mi ranging from 
(11.6, 13.2, 13.1) mJy to (25.7, 27.1, 35.8) mJy, depend- 
ing on the depth of the observation, with the deeper ob- 
servations being confusion limited. All of SERVS is cov- 
ered to Level 6 or deeper, with significant areas as deep 
as Level 3 (5cr limiting flux density «7mJy at 160 /im 
with the PACS 15 instru ment). Full details of the Her- 



MES survey are given in Oliver et al. (2012 1. 

The wide area component ot the SC U 13 A- 2 Cosmol- 
ogy Survey (S2CLS) will cover the XMM-LSS, Lockman, 
EN1 and CDFS fields. The survey will be performed at 
850 /im to a root mean square (rms) noise of 0.7 mJy. 

5.4. Radio surveys 

The Australia Te lescope Large Area Survey (ATLAS; 
Norris et al. 2006) overlaps with much of the SERVS 
fields in ESI and CDFS. The ATLAS survey will have an 
rms of 10 /LtJy and a spatial resolution of ss 8" at 1.4 Ghz 
across both fields (Banfield et al. in prep) . A prelimi 



nary release of the survey data for C DFS (Norris et al. 
2006] and ESI ( [Middelberg et al |2008| has an rms noise 
ot 30 /i Jy. An image showing the overlap between the 
SERVS and A TLAS fields can be f ound in the accompa- 
nying paper by Norris et al.| (2011a). Subsequent ATLAS 
data releases will be published shortly by Hales et al. 
(2012) and Banfield et al. (2012), both in preparation. 

The EN1 and Lockman fields have been surveyed at 
610 MHz with the G iant Meter Wave Telescope (GMRT; 
Garn et al. 2008a b, respectively). These surveys have a 
mean rms of « 60 /xJy at a spatial resolution of ss 6". 

The Very Large Array (VLA) has conducted several 
surveys in the SERVS fields. The Lockman field includes 
the deepest radio survey at 1.4 GHz to date, a single 
40' x 40' pointing centered at 161. 5d, +59.017d, reaching 
a 5 a detection limit of 15 /iJy near the center of the pri- 
mary beam (Owen fc Morri s on 2008), and further deep 



coverage by Ibar et al. ( 2009 ) . Simpson et al. ( 2006 ) have 



surveyed the Subaru Extragaiactic Deep Survey region 
with the VL A to a detection limit of 100 /iJy, and Bondi 
eT^(|2007] have surveyed the VVDS field in XMM-LSS 



aTT5oth~rJTDMHz with the GMRT and 1.4 GHz with the 
VLA to limits of « 200 and 80 /iJy, respectively. The 
Faint Images of t he Radio Sky at Twenty cm (FIRST. 
White et al.|[l997] survey covers the Lockman, EN1 sur- 



veys and part ot the XMM-LSS survey at 1.4 GHz to a 
sensitivity limit of ~ 1 mJy at a spatial resolution of « 5". 

Two pathfinder telescopes for the proposed Square 
Kilometer Array are currently under construction in the 
southern hemisphere. Both of these telescopes will un- 
dertake continuum surveys that will cover the south- 
ern SERVS fields. The Evolutio nary Map of the Uni- 



verse (EMU, Norris et al. 2011b) survey with the Aus- 
tralian Square Kilometer Array Pathfinder (ASKAP) 
will cover the whole southern sky to 10/iJy rms sen- 
sitivity at 1.4 GHz with a 10" FWHM (full width at 
half-maximum) synthesized beam. The South Af rican 
Karoo Array Telescope (MeerKAT, |Jonas| [2009| will 
conduct the MeerKAT International Giga-Hertz Tiered 
Extragaiactic Exploration (MIGHTEE) survey, which 
has a strong SERVS participation. MIGHTEE has 
several tiers, one of which will include the southern 
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Fig. 18. — Approximate 5<r depths of SERVS and near-IR and 
optical surveys covering the same areas. Ovcrplottc d is the SEP 
of a lGyr old stellar population at z = 2 from |Maraston| ( |2005| |. 
SERVS is in red, VIDEO in green, DXS in cyan, CFHT ^-bancTin 
orange and SpARCS in blue. Our target depth for optical bands 
imaging is shown as the magenta bar (several fields already have 
imaging to at least this depth). Note that no ground-based survey 
covers the entirety of the SERVS data, hence multiple surveys are 
overlaid. 



SERVS/VIDEO fields to 1 /iJy rms at 1.4 GHz with 
an w5" FWHM beam In the north the Low Fre- 
quency Array (LOFAR, Rottgering et al.||2011[ ) will tar- 
get the SERVS/SWIRE fields tor deep surveys at fre- 
quencies of ~ 100 MHz. In the northern hemisphere, 
the WODAN (Westerbork Observations of the Deep 
APERTIF northern-Sky) is being proposed for Wester- 
bork+APE RTIF and will match EMU sensitivity and 
resolution ( Rottgering et al.||2011 ). 



5.5. X-ray surveys 

T he XMM-LSS field o verlaps with the XMM-LSS sur- 
vey ( |Pierre et al.||2007D and the Subaru / XMM-Newton 
Deep' Survey (SXDS, [(Jeda et al.|[2008| . | Wilkes et alT 
( 2009 ) have a deep C handra survey overlapping with the 
deep VLA pointing of Owen & Morrison (2008J) in Lock- 
man. In EN1, Chandra program 690062 (P.l. Nandra) 
covers « 1 deg 2 of the central portion of the field. 

5.6. Further ground-based data taken or to be taken by 
the SERVS team 

Further multiwavelength ancillary data on the SERVS 
fields are currently being obtained with the overall goal of 
matching the SERVS depth in shorter wavebands. These 
data will be made available as part of the overall SERVS 
public release. Observations are concentrated on longer 
optical and near-IR wavebands as these are generally 



redshifts (van Dokkum et al. 2006 


Brammer et al.||2008 


Ilbert et al. 12009 Cardamone et al. 

in ii i m im — n-n ' i — • — i 1 ' 


2010j). In the optical, 



rower bands allow higher-fidelity photometri c re dshifts 



than the Johnson-Cousins system (see Figure 19 1 

Target depth in the optical is an AB magnitude « 25 
and ss 23 in the near-infrared (the area covered by_ 



VIDEO will be significantly deeper than this). Figure 18 
shows a comparison of the SERVS depth with those of 
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the other major optical/near-IR surveys planned or in 
progress. 

Observations with SuprimeCam on the Subaru tele- 
scope have been carried out in i and z bands as part of 
the Gemini-Subaru (PI A. Verma) and Keck-Subaru (P.I. 
S.A. Stanford) time swaps. The z-band observations are 
concentrated in the northern fields, as VIDEO will cover 
the southern fields in Z. 

An optical imaging campaign of the ESI and CDFS 
fields using the CTIO-4m Mosaic camera has been com- 
pleted (in November 2009 & October 2010, with a total of 
9 nights). The specific goals were (1) to complete imag- 
ing in r and i to 24.2 and 23.2 (Vega) over the whole area 
of the ESI and CDFS fields and (2) to obtain deeper r- 
and z-band data (to 25.0 and 24.0) in the center of ESI to 
complement deep surveys by VISTA and warm Spitzer, 
and where deep data at shorter wavelengths exist from 
ESO. Images are currently being processed. 

An ii-band imaging campaign, led by M. Lehnert has 
obtained ff-band data with WIRCAM on the CFHT, to 
match the UKIDSS DXS J and K imaging in the ENl 
field. These data are currently being analyzed. 

5.7. Photometric redshifts 

The wealth of ground-based data available in the 
SERVS fields in conjunction with the IRAC data ob- 
tained through the SERVS project provides the ideal ba- 
sis for a robust determination of photometric redshifts, 
stellar masses and other stellar population properties 
from SEDs. Here the photometric redshifts were com- 
puted via S EP fitting using the HyperZ code ( Bolzonella 
et al. 2000) and the |Maraston (2005) stellar popula- 
tion templates, following the procedure outlined in Pforr, 
Maraston & Tonini (2012b, in prep). In brief, the best 
fit is determined by minimizing the x 2 between a large 
grid of population model templates (for various star for- 
mation histories, ages, metallicities, reddening, and red- 
shifts) and the photometric data. Specifically, we use 
a template set with exponentially declining star forma- 
tion rates, metallicities ranging fr om 1/5 to twice solar , 
a minimum age of 0.1 Gyr and a Calzetti et al. (2000) 
reddening law with extinction parameter Ay ranging be- 
tween and 3 mag. This template set gives the minimum 
variance between photometric redshifts and the spectro- 
scopic redshifts calibration set. 

By comparing to available spectroscopic redshifts us- 
ing the compilation of redshifts in the SWIRE fields by 
Vaccari et al. (2012, in preparation), a redshift accuracy 
of Az/(1 + z) = 0.011 ± 0.072 is obtained for a dataset 
with an optimal wavelength coverage (i.e. U, g, r, i, Z, 
J, K, IRAC1 [3.6], IRAC2 [4.5] for the ENl field). Due 
to the small 16 number of SERVS objects with available 
spectroscopic redshifts, we further confirmed the redshift 
accuracy by cross-matching t he VVDS survey in th e ex- 
tended Chandra Deep Field (|Le Fevre et al.||2005|) with 



1600 n 



the IRAC SIMPLE survey rtDamen et al. 2009 [""obtain 



ing 831 matches with redshifts and 4.5 /im fluxes > 3/iJy. 
However, the diversity inherent to the ancillary data cov- 
erage for the SERVS fields impairs the redshift deter- 

16 In the ENl field, among the 2997 extended sources with opti- 
mal wavelength coverage (respectively 8513 with nonoptimal wave- 
length coverage) , only 30 (respectively 60) have spectroscopic red- 
shifts. 
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Fig. 19. — Photometric redshift distribution of SERVS galaxies 
in ENl (black histogram) . Photometric red shifts are o btained from 
SED fitting using the HyperZ code and the Maraston ( 2005 ) stellar 
population synthesis templates. Included are only objects classified 
as extended in the optical and near-IR bands in order to exclude 
stars and AGN-dominated sources. Additionally the red histogram 
shows the subsample of sources which are detected in all filter bands 
(U,g,r,i,Z,J,K, IRAC1, IRAC2 for ENl) which provides the most 
robust photometric redshift estimates. From Pforr, Maraston & 
Tonini (2012b, in preparation). 

mination particularly for objects with a narrower wave- 
length coverage. This is explored in detail in Pforr et 
al. (2012b, in preparation). In summary, photometric 
redshifts display larger scatter when the wavelength cov- 
erage used for the fitting is narrow, i.e. does not incorpo- 
rate the rest-frame UV, optical or near-IR rest-frame or 
a combination thereof, particularly when spectral breaks 
such as the 4000A break and the Lyman break are not 



2009 and 



covered by the filter setup (see Ilbert et al^ 
|Bolzonella et al.|2000| for stud ies on photometric redshift 
accuracies). 

The distribution of photometric redshifts obtained for 
extended 17 objects, i.e. galaxies, in ENl is shown in Fig- 
ure |19| The red histogram highlights the most robust 
sample (so-called "gold" sample) with the best wave- 
length coverage, i.e. corresponding to the case where 
the object is detected in all filter bands. The black his- 
togram includes also objects that were detected in less 
photometric bands. 



6. ONGOING SCIENCE WITH SERVS 

In this Section, we summarize already published or on- 
going work by SERVS team members to underline some 
of the key science goals of the SERVS survey. 



17 The decision whether an object is extended is based on stel- 
larity flags provided with the ancillary data in optical and near-IR 
filter bands. Sources with a pointlike nature indicative of stars or 
AGNs were excluded since the model template we use in the fitting 
does not include nonthermal emission 
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6.1. Galaxies and their environments 

Sampling ~ 0.8 Gpc 3 , SERVS is large enough to con- 
tain significant numbers of objects while still being deep 
enough to find L* galaxies out to z « 4 (Figure 16). 



By combining the five different fields of SERVS, the sur- 
vey effectively averages over large-scale structure, and is 
able to present a true picture of the average properties 
of galaxies in the high redshift Universe. 

Galaxy-galaxy correlations are being computed by van 
Kampen et al. (2012b, in preparation) in the SERVS 
fields. The five large, well-separated, SERVS fields en- 
able us to average out the effects of large-scale struc- 
ture on such measurements. Initial results on the EN1 
field show the evolution of the correlation function be- 
tween high redshifts (z > 1.3) and intermediate redshifts 
(~ 0.8) using simple [3.6] — [4.5] color cuts. 

SERVS is deep enough and wide enough to find field 
clusters at z > 2, should they exist in significant num- 
bers. The SERVS fields lie within the SpARCS fields 
(see survey description in Section 5.1 DeGroot et al. 
2012, in prep), The deep SERVS observations will al- 
low both a more accurate measurement of the faint end 
of the luminosity function of known SpARCS clusters 
which fall within the SERVS footprint, and also the de- 
tection of new clusters at higher redshifts than possible 
with SpARCS. 

In addition, Geach et al. (2012, in preparation) are 
pursuing a cluster selection technique using photomet- 
ric redshifts combined with Voronoi tessellation in an 
attempt to identify further, mostly lower mass, cluster 
candidates. 

6.2. The Active Galactic Nuclei engine 

The [3.6] and [4.5] bands are important diagnostics of 
AGN SEDs, as they are where host galaxy light and hot 
dust emission from the torus overlap in the SEDs of many 
dust-obscured AGN and quasars at moderate redshifts 
(z ~ 1). In unobscured, or lightly obscured objects, 
this is where the optical/UV emission from the ac c retion 
disk transitions to the hot dust emission. Petric (20101 
& Petric et al. (2012, in preparation) present SEDs of 
AGNs and quasars selected in the mid-infrared, and use 
SERVS data to help apportion the different sources of 
near-infrared light. The luminosities of the hosts them- 
selves, if free from contamination by AGN-related light, 
can also be used to study the stellar masses of the host 
galaxies. 

6.3. AGN and their environments 

Current models for galaxy formation indicate that 
AGN and quasar acti vity play an import ant role in 
galaxy formation (e.g. |Hopkins et al. 20061, regulating 
the growth of their host galaxies through feedback (see 
for example |Schawinski et al. 12007} or more recently |Far-| 
|rah et al.|2012| who provide new and strong evidence for 
AGN feedback). However, the exact nature of this feed- 
back process is unclear. Environments in which AGN 
and quasars lie can indicate the masses of the dark halos 
they inhabit, and also ho w these masses dep end on AGN 
luminosity and redshift (Farrah et al. 2004). These can 
illuminate models for feedback, tor example, a prepon- 
derance of AGN in massive halos, accreting at relatively 
low rates might be an indicator that their host galaxies 
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Fig. 20. — Stacked source over-density vs radial distance for the 
11 QSOs in the redshift range of 2.8 < z < 3.8. The first bin 
has a radius of 700 kpc and the other bins are of the same area 
as the first. The error bars show the Poisson error on the number 
counts. The dashed line shows the subtracted local background 
level (zero level) determined from an annulus of 2Mpc 2 (400") 
from the QSOs. The dotted line shows, for comparison, the global 
background as determined from taking the average source density 
in large apertures over the SERVS fields. This is the s ource density 
before being corrected for completeness. Figure from Faldcr et al. 
(2011) ■ 



are no longer growing rapidly ( |Hopkins et al.|[2007 ). At 
low redshifts (z < 0.6) the SPSS h as been used to suc- 
cessfu lly perform these experiments ( Padmanabhan et al. 
2009[ ). SERVS is able to take these studies to z » 1. 

One particular area where SERVS is uniquely valuable 
is in determining th e environments of high redshift AGN. 
Falder et al. (2011 ) find significant (> 4er) overdensity of 
gaiaxies around QSOs in a redshift bin centered on z~ 
2.0 and an (> 2 a) overdensity of galaxies around QSOs 
in a redshift bin centered onz~ 3.3 (see Figure 20). 

Nielsen et al. (2012, in preparation) are investigating 
the environments of AGN and quasars selected in the 
mid-infrared. For the first time the environments of lu- 
minous quasars at 0.8 < z < 3 are being characterized, 
enabling a comparison of the environments of dust ob- 
scured and normal quasars at these redshifts. 

6.4. High-z quasar searches 

The unique multiband SERVS data set will be a valu- 
able tool for constraining the faint end of the quasar lumi- 
nosity function at high redshifts. Quasar searches have 
been or will be carried out in the SERVS fields using 
the combination of SERVS, DXS, VIDEO and SpARCS 
data, in addition to the SERVS CTIO and Subaru an- 
cillary data. This large range of wavelengths allows 
for the rejection of many contaminants of high-z quasar 
searches on the basis of near-infrared photometry alone. 
Current es timates of the high-z quasar luminosity func- 
tion (e.g. Willott et al. 2010 suggests somewhere be- 
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Fig. 21. — Two representative infrared-faint radio sources 
(IFRS). The greyscale is the [3.6] SERVS data, and the contours are 
the 20 cm image, with contour levels of (1, 2, 3, 4, 5) mjy/beam. 
The left hand image is a non-detection(CS0194) and t he right-hand 
image i s a candidate detection (CS0114), Figure 2. of |Norris et al.| 
1201 la) ■ 



tween 3-14 quasars at 5.5 < z < 6.5 and 1-5 quasars at 
6.5 < z < 7.5. The large range is due to the uncertainty 
in the faint-end slope of the quasar luminosity function, 
and SERVS will be able to constrain this well. 



6.5. Infrared-faint radio sources 
Norris et al. (120061) describe the discovery of infrared- 



faint radio sources (IFRS), a population of radio sources 
with hos t galaxy fluxes well be low the limit of the SWIRE 
survey. Huynh et al. ( 2010| ) used deep IRAC data to 
place even more stringent limits on them, and concluded 
that these are most likely radio-loud AGN wit h faint 
host galaxies, but the sample to date is small. |Norris| 
et al. (2011a) present an initial study of this hitherto 
unsuspected population with SERVS, including stack- 
ing of objects that are too faint to be detected, even in 
SERVS, and that may represent very high redshift radio- 
loud galaxies, possi bly s uffering from significant dust ex- 
tinction (see Figure 21). 

Using a combination of SERVS data and GMRT/VLA 
radio observations of the Lockman Hole at 610 MHz and 



1.4 GHz, Afonso et al. (2011) study a sample of Ultra 
Steep Spectrum (USS) radio sources and suggest the 
likely existence of higher rcdshifts among the submilli- 
jansky USS population, raising the possibility that the 
high efficiency of the USS technique for the selection of 
high-redshift sources remains even at the sub-mJy level. 



6.6. Obscured star formation 

Although SERVS cannot be used as a direct indicator 
of obscured star formation on its own, overlap with sur- 
veys by SCUBA-2 (S2CLS) and Herschel (HerMES) will 
allow more reliable source identification than possible us- 
ing shorter wavelength data, and better characterization 
of any extinction of the stellar light, as well as the stel- 
lar mass of the galaxy. Based on recent simulations of 
Herschel observations, we expect to detect ~ 700 uncon- 



fused sources per d eg in > 3 Herschel bands (Fernandez> 



Conde et al. 2008). At the SERVS depths, we expect to 



detect > 95% of these sources in both IRAC bands, and 
with the aid of our ancillary optical and near-infrared 
data, obtain photometric redshifts and stellar masses for 



~ 12000 sources. This will be sufficient to study trends 
in star formation rate with stellar mass and redshift, for 
example, to test the idea of "downsizing" of the most 
actively star-forming galaxies. 



7. SUMMARY 

The Spitzer Extragalactic Representative Survey 
(SERVS) is designed to open up a medium-depth, 
medium-area part of parameter space in the near- 
infrared, covering 18 deg 2 to «2yuJy in the Spitzer [3.6] 
and [4.5] bands in five highly observed astronomical fields 
(EN1, ESI, Lockman Hole, CDFS and XMM-LSS). The 
five SERVS fields are centered on or close to those of 
corresponding fields surveyed by the shallower SWIRE 
fields, and they overlap with several other major sur- 
veys covering wavelengths from the X-ray to the ra- 
dio. Of particular importance are near-infrared data, 
as these allow accurate photometric redshifts to be ob- 
tained. SERVS overlaps exactly with the 12 deg 2 of the 
VISTA VIDEO survey in the south, and is covered by 
the UKIDSS DXS survey in the north. SERVS also has 
good overlap with HerMES in the far-infrared, which cov- 
ers SWIRE and other Spitzer survey fields, with deeper 
subfields within many of the SERVS fields. Sampling 
~ 0.8 Gpc 3 and redshifts from 1 to 5, the survey is large 
enough to contain significant numbers of rare objects, 
such as luminous quasars, ultra- luminous infrared galax- 
ies, radio galaxies and galaxy clusters, while still being 
deep enough to find L* galaxies out to z w 4. In this arti- 
cle, we have described the Spitzer observations, the data 
processing and the wealth of ancillary data available in 
the fields covered by SERVS. Mosaics and catalogs will 
be made available to the community in the summer of 
2012 through the Infrared Science Archive (IRSA). 
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APPENDIX 



NUMBER COUNTS FOR THE EN1. ESI, LOCKMAN & CDFS FIELDS 



Number counts were derived for the five fields. Source colors were also compared to 2MASS (see details in Section 4.2 ). 
We present here the four remaining fields: EN1, ESI, Lockman and CDFS. All of them show consistent number count 
features and completeness levels. SERVS bright stars colors for all of the fields are coherent with 2MASS. 





Fig. A.l. — Left side: SERVS number counts versus flux at [3.6] and [4.5] for the EN1 (top) and ESI (bottom) fields. Grey dashed 
lines show the selection of sources used in the bottom plot (with fluxes as 0.3 < / < 2mjy). Right side: Color-color plot showing 
^"2MASS — A^3.6jtm(Vega) versus M^.s^m (Vega) — M4.5 Mm (Vega) for sources within the flux range defined by the grey dashed lines above, 
plus a cut in stellarity index > 0.95 and the existence of a 2MASS K-band measurement as an additional constraint. Red dashed lines help 
pinpoint the location of the (0,0) point in this diagram. 
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Fig. A. 2. — Left side: SERVS number counts versus flux at [3.6] and [4.5] for the Lockman (top) and CDFS (bottom) fields. Grey 
dashed lines show the selection of sources used in the bottom plot (with fluxes as 0.3 < / < 2mjy). Right side: Color-color plot showing 
^"2MASS — A^3.6M"i(Vega) versus M^.s^m (Vega) — M4.5 Mm (Vega) for sources within the flux range defined by the grey dashed lines above, 
plus a cut in stellarity index > 0.95 and the existence of a 2MASS K-band measurement as an additional constraint. Red dashed lines help 
pinpoint the location of the (0,0) point in this diagram. 



SEXTRACTOR & MOPEX PARAMETERS USED FOR THE IMAGE PROCESSING AND THE CATALOGS EXTRACTION 



Mopex was u sed t o process and co-add the raw images and SEXTRACTOR to extract information on the sources 
within. In Table [O] we list the core parameters used in both. 
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TABLE A.l 

Values of the more important Mopex & SExtractor parameters 



Program (module) 


Parameter 


Value 


t\ T/"itz>ity /TMT r ri7r^ r r 'i 

SyXKJrhijL LZi 1 HiV^ 1 ) 


Detection_iVlax_Arca 


1 on 


IVIOPEX ^Uri/ 1 fjKj 1 J 


Detect ion _Min _Arc a 


U 


MOPEX (JJrjlri/ClJ 


Detect ion_Threshold 


i 

1 


MOPEX (MObAlGliN 1 J 


TMTTPTD T) A T "C 1 T XT r\ TA 

UN 1 rjKF.Mrj 1 HOD 


1 


MOPEX (MObAlGDUALOU ILil^K) 


TV TT\T /~\ T Trpr T1 t a /"ITT! 


2 


TV J i-\-r-iT-iir /ATAO A T/TAT7 A T /"A T TT^T TTTTD \ 

MOPEX (MObAlODU ALOU ILlrLKj 


MIJNLOU IL-t 4 HAO 


0.51 


MOPEX (MUbAlGUU ILliiiK) 


r PTJ"D T7 O TJT rMDTT AM 

1 HKJ^bri_(J.r 1 1(J1N 


1 


MOPEX (MObAlGOU ILlhiK) 


BUI l(JM_lriKrjbri(JLD 





MOPEX (MUoAlbU U 1 LlJiiK J 


T"^rAT3 r PXJ'D TTQTJT^T TA 





MOPEX (MObAlGOU iLlh/K) 


A JTAT TDTA7" TVTT TA /T 


3 


MOPEX (MUoAlOHMAoK) 


A TTTVT f~*r~\\ TT?T> A T? 


4 


A /TrtnTiv /'AyT/'AO A T /~i "D A T A CIV'S 

MOPEX (MUoAlOHMAoK) 


jj A v rini 7'T7 , T3 A T7 1 

M A X _0 (J V ri/ H A O ii/ 


100 


SExtractor 


TAT7 ir rrr , ^ r r attm a "dtt 1 a 
JJili 1 ri/L/ 1 _M1IN AKri/A 


o . U 


O-CjA 1 KAt 1 (JK. 


UHi X TZj\_' J L JT1 i\,J_jOrl 


U.l 


SExtractor 


ANALYSIS_THRESH 


0.4 


SExtractor 


FILTER 


Default 


SExtractor 


DEBLEND_NTHRESH 


64 


SExtractor 


DEBLEND_MINCONT 


0.005 


SExtractor 


SEEING_FWHM 


2.0 


SExtractor 


BACK_SIZE 


32 


SExtractor 


BACK_FILTERSIZE 


5 


SExtractor 


BACKPHOTO_TYPE 


LOCAL 


SExtractor. 


WEIGHT_TYPE 


MAPJWEIGHT * 



Notes: Similar values can be found in |Lacy et eU7| ( |2005[ ) and |Ashby et al.| | |2009| | 

* Using coverage map as inverse variance weight map 
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